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^ . ABSTRACT 

Context. Since its launch, the Fermi satellite has firmly identified 5 pulsar wind nebulae plus a large number of candidates, all powered 
by young and energetic pulsars. HESS J1857+026 is a spatially extended y-ray source detected by H.E.S.S. and classified as a possible 
pulsar wind nebula candidate powered by PSR J1856+0245. 

Aims. We search for y-ray pulsations from PSR J1856+0245 and explore the characteristics of its associated pulsar wind nebula. 
Methods. Using a rotational ephemeris obtained from the Lovell telescope at Jodrell Bank Observatory at 1 .5 GHz, we phase-fold 
36 months of y-ray data acquired by the Large Area Telescope (LAT) aboard Fermi. We also perform a complete y-ray spectral and 
C*"} ■ morphological analysis. 

Results. No y-ray pulsations were detected from PSR J1856+0245. However, significant emission is detected at a position coincident 
\Q ' with the TeV source HESS J1857+026. The y-ray spectrum is well described by a simple power-law with a spectral index of T = 

1.53±0.11 st at±0.55 syst and an energy flux of G(0. 1-100 GeV) = (2.71 ±0.52 st!lt ± 1.51 syst )x 10"" ergs cirT 2 s -1 . The y-ray luminosity 

is L y PWN (0.1-100 GeV) = (2.5 ± 0.5.„„, + 1.5 IV ,„) x 10 35 (5^) ergs s assuming a distance of 9 kpc. This implies a y-ray efficiency 
of ~ 5 % for E = 4.6 x 10 36 erg s' 1 , in the range expected for pulsar wind nebulae. Detailed multi-wavelength modeling provides new 
constraints on its pulsar wind nebula nature. 

Key words, pulsars : general, pulsars : individual object : PSR J1856+0245, ISM : individual object : HESS J1857+026, Gamma 
% , rays : general 

1. Introduction are ~ 1%, consistent with TeV observations (lAckermann et al.l 

SoTD)- 

Pulsar wind nebulae (PWNe) are bubbles of shocked relativistic 

particles produced by the interaction of the pu lsar's wind with The presence of a pulsar close to the source position is an 
the surrounding medium dGaensler el al.ll2006l) . Since 2003, the important clue to identify a PWN, which often requires infor- 
continuous observations of the Galactic Plane by Cerenkov tele- mation from the radio/X-ray wavelengths. Radio/X-ray PWNe 
scopes have yielded the detection of more than 60 Galactic TeV are often associated with TeV extended sources offset from their 
sources. Among them, PWNe are the dominant class with 29 pulsar s, which can be exp lained by an inhomogeneous environ- 
firm identifications. In the GeV energy range, 7 PWNe have ment dHintonetalJ 12010) . In such sources, TeV radiation can 
been firmly identified by the Fermi-LAT, all of them having their be explained by IC scattering of accelerated leptons on ambient 
Inverse Compton (IC) peak at energies higher than 100 GeV. photon fields (CMB, IR, ...) or by 7r° decay from the interaction 
They are all powered by energetic pulsars and their y-efficiencies of accelerated hadrons with nuclei of the interstellar medium. 
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HESS J 1857+026 is a very high energy (VHE) y-ray 
source detected by H.E .S.S. during the Galactic Plane Survey 
dAharonian et all 120081) . The extended (~ 0.11°) TeV source 
was identified as a PWN candidate after the discovery of PSR 
J 1856+0245 (offset ~ 0.12°) in the Arecibo PALFA survey 
dHessels et alj|2008l) with a Dispersion Measure (DM) of 622 
cm -3 pc. Recently, MAGIC reported a measured extension in 
the 0.2-1 TeV energy range significantly larger (0.22°) than the 
extension reported b y H.E.S.S. in the 0.6-80 TeV energy range 
(lKlepseretalJl201ll) . PSR J 1856+0245 is an energetic pulsar 
(E = 4.6 x 10 36 erg s ) located in a crowded region, 1.3° from 
the bright SNR W44 (Abdo et al. 2010) and 0.6° from the fainter 
SNR HESS J1858+020 o n which only an up per limit could be 
set using Fermi-LAT data dTorres et al.ll2.Ql lb . Significant emis- 
sion coincident with HESS J 1857 +026 was observed abov e 100 
GeV using Fermi-LAT observations dNeronov et al.ll20lol) . 

Here, we report in detail GeV observations of the HESS 
J1857+026/PSR J1856+0245 system using Fermi-LAT observa- 
tions and discuss their implications for the nature of the source. 



2. LAT description and data selection 

The LAT is a y-ray telescope that detects photons by conversion 
into electron-positron pairs and operates in the energy range be- 
tween 20 MeV and 300 GeV Details of the instrument and data 
processing are given i n lAtwood et al.l d2009l) . The on-orbit cali- 
bration is described in I Abdo et al.l d2009al) . 

The following analysis was performed using 36 months of 
data collected from August 4, 2008 to August 31, 201 1 within a 
10 x 10° square around the position of HESS J1857+026 aligned 
with Galactic coordinates. We excluded y-ray s coming from a 
zenith angle larger than 100° because of possible con tamination 
from secondary y-rays from the Earth's atmosphere (Abdo ^etaD 
2009b). We used the P7 V6 Instrument Response Functions 
(IRFs), and selected the 'Source' events which correspond to the 
best compromise between the number of selected photons and 
the charged particle residual background for the study of point- 
like or slightly extended sources. 



3. Data analysis 

3.1. Timing analysis of PSR J 1856+0245 

With its large spin-down power, PSR J 1856 +0245 is one of the 
more energetic radio pulsars known. Its spin period of 80.9 ms 
and characteristic age of 20.6 kyr are similar to those of the 
Vela pulsar. The DM and NE2001 electron density model of 
the Galaxy assign PSR J 1856+0245 a distance of ~ 9 kpc 
dCordes et al.ll2002l) . 

This pu lsar is not monito red as part of the LAT pulsar timing 
campaign ( Sm ith et all2008l) . as it was discovered subsequently, 
but has nevertheless been regularly o bserved with the L ovell 
telescope at Jodrell Bank Observatory (iHobbs et al.l l2~004). The 
ephemeris of PSR J 1856 +0245 used in the analysis was obtained 
using 82 observations at 1.5 GHz made with the Lovell tele- 
scope between May 4, 2008 and September 4, 201 1. The event 
arrival times were corrected to the Solar System Barycenter us- 
ing the JPL DE405 Solar Syste m ephemeris. The TEMP02 tim- 
ing package (Hobbs et al. 2006) was then used to build the tim- 
ing solution. We fit the radio times of arrival (TOAs) to the pul- 
sar rotation frequency and first four derivatives (in order to re- 
move timing noise). By including a fourth derivative the RMS 
of the timing residuals decreases by about 30%, from 1.55 ms 
to 1.08 ms, equivalent to a decrease from 19 to 13 milliperiods. 
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Period epoch (MJD) 


55128 


Start time (MJD) 


54615 


End time (MJD) 


54570 


Number of TOAs 


82 


TOA rms (ms) 


1.08 



Table 1. Parameters of the fit of the TOAs. v correspond to 
the rotational frequency of the pulsar and v w its time derivative 
of order i. In parentheses are the 1 <x uncertainty on the least- 
significant digits quoted, from TEMP02. 



The reduced^ 2 of the fit also decreases significantly, from 22 to 
5. The fit results are summarized in Table [TJ This timing solu- 
tion will be made available through the Fermi Science Support 
Cente£](FSSC). 

For the LAT analysis, photons with energies above 100 
MeV and within a radius of 1.0° of the radio pulsar posi- 
tion a(J2000) = 18 /! 56 m 50.937 ! , c5(J2000) = +02°45' 47.046" 
were selected using an energy-dependent cone of radius 6 < 
max(5.12° x (£7100 MeV) ,0.2°) and phase-folded using the 
radio ephemeris previously described. This choice takes into 
account the instrument performance and improves the single- 
photon signal-to-noise ratio over a broad energy range. No sig- 
nificant pulsation was detected for all tested energy bands (100 
MeV - 300 GeV, 100 MeV - 300 MeV, 300 MeV - 1 GeV, > 
1 GeV). Following the procedure used by Roma ni et al 1 (l2011h . 
we fitted a point source at the position of PSR J 1856+0245 in 
the 0.1-1 GeV energy range assuming a power-law of index 
1.62 and a cut-off energy at 2.8 GeV to derive a 99% Bayesian 
upper limit on the flux of 3.27 x 10~ 8 ph ctrT 2 s _1 , well below 
typical y-ray fluxes reported for pulsars detected by Fermi-LAT 
(Abdo etaUl2010ch . 



3.2. Spatial and spectral analysis 

Two different tools were used to perform the spatial and spec- 
tral analysis: gtlike and p ointlike, gtlike is a binned 
maximum-likelihood method (M attox et al.ll 19961) implemented 
in the Science Tools distributed by the FSSC. pointlike is 
an alternate binned likelihood technique, optimized for char- 
acterizing the extension of a sourc e (unlike gt like), that was 
extensively tested against gtlike dKerrll201lb . These tools fit 
a source model to the data along with models for the instru- 
mental, extragalactic and Galactic components of the back- 
ground. In the following analysis, the Galactic diffuse emission 
is modeled by the standard LAT diffuse emission ring-hybrid 
model galJ2yearp7v6_v0.fits. The residual cosmic -ray back- 
ground and extragalactic radiation are described by a single 
isotropic component with a spectral shape described by the 
file iso_p7v6source.txt. The models have been released and de- 
scribed by the Fermi-LAT Collaboration through the FSSQ^. 

The source significance is measured by a test statistic (TS) 
defined as TS= 2 (log (L\ ) - log {L$)), where L\ corresponds to 



FSSC: http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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the likelihood obtained by fitting a model of the source of inter- 
est and the background model and Lo corresponds to the likeli- 
hood obtained by fitting the background model only. In the fol- 
lowing, the correspondence between the significance and the TS 
value is evaluated from the^ 2 distribution with 4 degrees of free- 
dom (position and spectral parameters). 

The 41 sources within 15° of HESS J1 857+026 in the 
Second Fermi-LAH catalog dNolan et alj|2012l) were taken into 
account. We refitted all spectral parameters of the 16 sources 
within 5° around HESS J1857+026. The region includes the 
bright SNR W44, known to interact with its environment. 
Extended and only 1.3° from HESS J1857+026, W44 could in- 
fluence our fit. We refitted it assuming a n elliptical r ing an d 
obtained results consistent with those of Abdoetal. (2010a). 
The centroid is located at a = 18 ft 56 m , 6 = +01°22'. The fit- 
ted semi-major and semi-minor axes are respectively maj /2 = 
(0.33 + 0.l0 stat )°, min/2 = (0.20 + 0.02 sfflf )°. The angle of the 
semi-major axis from celestial North, taken positive toward in- 
creasing Right Ascension, is (327 + 22 stat y. 

3.2.1 . Shape and position of HESS J1857+026 counterpart 

Source shape analysis requires the best possible angular resolu- 
tion. Since the source has a hard spectrum (see Section l3.2.2l ) we 
made a compromise between statistics and resolution by select- 
ing photons above 10 GeV. This drastically reduces the contribu- 
tion of the Galactic diffuse background and improves the single- 
photon angular resolution. Fig. 1 (Top) presents a LAT TS map 
in the energy range of 10 GeV to 300 GeV. To each pixel is asso- 
ciated a TS value calculated assuming a point source in its cen- 
ter and fitting only the flux of the source assuming a power-law 
spectrum with a spectral index of 2. A source coincident with 
HESS J1857+026 is clearly visible. We determined the exten- 
sion of the source using pointlike with three different models 
: a point source, a uniform disk and a Gaussian. No significant 
extension was obtained above 10 GeV. The GeV emission was 
fit to position a(72000) = 18 /2 57" ! , 5(72000) = +02°45' with 
an average statistical error of 0.05°, consistent with the H.E.S.S. 
position, a(72000) = 18 A 56 ffl 50.80\ 5(72000) = +02°45'50.2". 

3.2.2. Spectral analysis 

Fig. 1 (Middle) shows a TS map of the region in the energy range 
0.1-1.3 GeV. There is excess emission near HESS J1857+026 
located at a(72000) = 18 /, 54" ! , 5(72000) = +02°59'. This ex- 
cess is inconsistent with that of HESS J1857+026 and was added 
to the background model. This additional background source 
was fitted assuming a pure power-law with an integrated flux 
of F(0. 1-100 GeV)=(2.31 ±0.37 stat ) x 10~ 7 ph cnr 2 s"\ a spec- 
tral index of F = 3.18 ±0.56 st at, which gives a significance above 
300 MeV of ~ 3.6cr (TS=20). 

Spectral analysis was performed using gtlike, selecting 
only 0.3-300 GeV to avoid the low energy range that is dom- 
inated by the diffuse Galactic background and subject to large 
systematic uncertainties. In this energy range, HESS J1857+026 
is well described by a pure power-law with an integrated flux 
extrapolated down to 100 MeV of F(0. 1-100 GeV)=(5.78 + 
0.85 stat ± 3.1 l sy st) x 10~ 9 ph ctrT 2 s _1 , a spectral index of 
r = 1.53 ± 0.11 stat ± 0.55 syst and an energy flux of G(0.1-100 
GeV)=(2.71 ± 0.52 stat ± 1.51 syst ) x 10" 11 ergs cnr 2 s"\ which 
gives a significance above 300 MeV of ~ 5.4cr (TS=39). The 
residual TS map after fit in Fig.l (Bottom) shows no significant 
excess. 



Fenni-LAT spectral points for HESS J1857+026 were ob- 
tained by splitting the 0.3-100 GeV range into 4 logarithmically- 
spaced energy bins plus a bin between 100 and 300 GeV which 
contains 22 photons corresponding to a TS of 14, shown in Fig. 
2. A 99 % C.L. upper limit is computed when TS<10 using the 
approach of N olan et al.l(l2012l) . The errors on the spectral points 
represent the statistical and systematic uncertainties added in 
quadrature. 

Three main systematic uncertainties can affect the LAT flux 
estimate for a point source: uncertainties in the Galactic diffuse 
background, uncertainties on the effective area and uncertainties 
on the source shape. The dominant uncertainty at low energy 
comes from the Galactic diffuse emission, estimated by chang- 
ing the normalization of the Galactic diffuse model artificially 
by +6% as in (Abdo 2010 b). Since it is computed for P6 IRFs, 
the 6% factor overestimates the bias. The second systematic is 
estimated by using modified IRFs. The fact that we do not know 
the true y-ray morphology introduces a last source of error. We 
derived an estimate of the uncertainty on the source shape by us- 
ing the best Gaussian model obtained by H.E.S.S. We combine 
these various errors in quadrature to obtain our best estimate of 
the total systematic error at each energy, which we propagate 
through to the fit model parameters. 

Assuming a distance of 9 kpc and isotropic emission, the y- 
ray flux corresponds to L y pwN (Q. 1-100 GeV) = (2.5 ± 0.5.5, n , ± 

1.5 JVJf ) x 10 35 (5^) ergs s _1 . Using the pulsar's E, yields a y- 
ray efficiency of ~ 5 %, one of the highest PWN efficiencies ob- 
served at GeV energies ( Ackermann et al. 2 01 lb . It is still in the 
range of expected values for PWNe seen by the F ermi-LAT and 
is clo se to the estimate of 3% using H.E.S.S. data (Mattana et al. 
I2009h . 

4. Supporting X-Ray measurement 

To find the flux of a potential X-ray PWN associated with PSR 
J 1856+0245, we analyzed a 39-ks Chandra ACIS-I observa- 
tion from February 28, 2011 (Obs. ID 12557). These data were 
recorded in the VFAINT and Timed Exposure (TE) modes and 
were analyzed using CIAC0 version 4.3.1 with CALDB 4.4.3. 
PSR J 1856+0245 is clearly detected as a point source, but there 
was no immediate evidence for extended emission surrounding 
this position. Based on a 30-ks XMM-Newton EPIC pn obser- 
vation (ObsID 0505920101), the unabsorbed flux of the pul- 
sar in the 2-10 keV range, assuming a power-law spectrum, is 
8.3^ 2 „ X 10~ 14 ergs s _1 . A more in-depth analysis of the X-ray 
properties o f PSR J1856+0245 as well as the XMM data will be 
presented in Bogda nov et al.l din pr eparation). 

Given that the size of the potential X-ray PWN is not known, 
we investigated an appropriately sized extraction region to see 
whether it contains statistically significant count excess above 
background. The extraction region is an annulus extending from 
2" - 15" from the position of the pulsar. The inner radius was 
chosen to avoid contamination from the pulsar. The outer ra- 
dius was chosen based on the X-ray PWNe observed for pul- 
sars with comparable E by scaling th eir angular size with their 
distance dKargaltsev & Pavlovll2008h . The background regions 
were chosen from several other source-free regions near the pul- 
sar. For the 2" - 15" extraction region we find an upper limit 
on the unabsorbed flux of 5 x 10~ 14 erg s _! cirT 2 (l-10keV, 3<x 
confidence), corresponding to a luminosity of 5 x 10 32 erg s _1 . 

3 C handra Interactive Analysis of Observations dFruscione et all 
120061) . 
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The counts in this region show a 2<x excess from zero counts. 
Assuming instead the TeV position and a 6' extraction region 
consistent with the H.E.S.S. morphology, we find an upper limit 
on the unabsorbed flux of 2 x 10~ 12 erg s _1 cirT 2 (1-lOkeV, 3<x 
confidence), corresponding to a luminosity of 2 x 10 34 erg s _1 . 
Given the marginal significance of the count excesses derived, 
we cannot convincingly claim the detection of a weak X-ray 
PWN. These luminosity limits are derived from the 3<x upper 
bound on the net count rate and assume a typical power-law 
spectrum of index 1 .5 for the PWN, a distance of 9 kpc, and a 
column density Nh = 4 x 10 22 crrT 2 based on the spectroscopic 
analysis of the XMM data. 

5. Discussion 

The nature of HESS J 1857+026 remains unclear. From the ob- 
servations presented here, we know that the GeV source is po- 
sitionally and spectrally consistent with the TeV source, sug- 
gesting a physical relationship. The limits on an X-ray PWN 
indicate a low magnetic field for any leptonic model, because 
a larger field would produce strong X-ray synchrotron emission. 
The spatial extent of the TeV source and the proximity to a pul- 
sar suggests a PWN. The distance, characteristic age, and energy 
available from the pulsar are known, although any distance es- 
timate based on DM has a significant uncertainty. Using these 
inputs, we investigate whether a plausible PWN model consis- 
tent with all the observations can be found. 

To investigate the global properties of the PWN, we apply 
a one-zone time dependent Spectral Energy Distribution (SED) 
model which r eproduces the mul ti-wavelength measurements 
from MAGIC dKlepser et al.ll2011l) . H.E.S.S. dAharonian et alj 
2008), as well as the LAT results and the 2" - 15 " X-ray upper 
limit d escri bed above. This m odel is described in lGrondin et alj 
(1201 lh and lAbdo et all d2010h . with a more detailed description 
(applied to a m ulti-zone scenario rather than one-zone) found in 
Section 3.10 of I Van Etten & Romanil d2011l) . 

The model computes SEDs from evolving electron popula- 
tions over the lifetime of the pulsar in a series of time steps with 
the energy content of the injected particle population varying 
with time according to the pulsar spin down. At each time step, 
new particles are injected and synchrotron, inverse Compton, 
and adiabatic cooling is computed. Synchrotron and inverse 
Compton fluxes are computed from the final electron spectrum, 
magnetic field, and photon fields. Starting assumptions include: 
the source of all particles and magnetic field in the nebula is 
PSR J 1856+0245, ambient photon fields are uniform and do not 
vary over the evolution time of the electron populations, and the 
nebula distance is 9 kpc. 

The large size of the nebula (~ 20 pc at 9 kpc) implies a 
middle-aged PWN, whose parent SNR has likely evolved from 
the free-expansion phase to the Sedov phase. The Sedov phase 
is expected to occur on a timescale of ts e dov ~ 3 kyr for an 
explosion of 10 51 erg, an ejecta mas s of 1OM0, and an am - 
bient medium density of 0.1cm 3 dTruelove & McKeeH l999). 
Eventually, the inward moving SNR reverse shock collides with 
the expanding PWN, which can happen as late as 5 time s the 
transition to the Sedov phase (Ivan der Swaluw et al.l l2004). The 
interaction of the PWN and the SNR reverse shock compresses 
the PWN, resulting in an increased magnetic field. Even in the 
spherically symmetric case the evolution is complex, with the 
PWN undergoing a series of oscillations due to the reverse shock 
interaction. The significant offset of the pulsar from the y-ray 
centroid implies either a PWN expansion into an inhomogenous 
medium, or an asymmetric reverse shock interaction. Either way, 



models such as Gelfand "et al.l d2009l) which compute the PWN 
compression assuming spherical symmetry cannot be applied. 

Given the lack of multi-wavelength data for HESS 
J 1857+026 and the complexity of the SNR-PWN interaction, 
we see little reason to adopt a complicated spatial model and 
therefore assume a simple spatial evolution. During the free- 
expansion phase of the PWN (assumed to be ~ 3 kyr) we adopt 
an expansion of R oc t, following which the radius evolves as R oc 
f 3 , appropriate for a P WN expanding in pressure eq uilibrium 
with a Sedov phase SNR dRevnolds & Chevarierll 19841) . We also 
adopt a simple magnetic field model, and assume that over the 
pulsar lifetime the magnetic field evolves as B oc t , following 
~ 500 years of constancy. This magnetic field evolution is sim- 
ilar to the r 1 3 - r 2 behavior d erived by Reynolds & Chevalierl 
d!984l) and the r 1 7 evolution of lGelfand et al.l d2009l) during the 
initial expansion phase. The power-law e volution of the mag - 
netic field is s imilar to that adopt e d by IZhang et al] (2008), 
iLemiereetafl d2009l) . iMaver et all (1201 2|) for similarly aged 
PWNe. 

Model fitting is achieved by minimizing the x 1 between 
m odel and da t a usin g the downhill simplex method described 
in iPress et al . (1992). For each ensemble of N variable param- 
eters we evolve the system over the pulsar lifetime and calcu- 
late^ 2 between model curves and flux data points. The simplex 
routine subsequently varies the parameters of interest to mini- 
mize the fit statistic. We estimate parameter errors by computing 
X 2 for a sampling of points near the best fit values and using 
these points to fit the A^-dimensional ellipsoid describing the 
surface of Ay 2 = 2.71. Assuming that the adopted PWN model 
is correct, and that errors are Gaussian, the projected size of this 
Ay 2 =2.71 ellipsoid onto each parameter axis defines the 90% 
multi-parameter (projected) error. 

Using the radio pulsar ephemeris, we tried to derive a value 
for the braking index using Eq.[T]where v, v w and n represent re- 
spectively the rotational frequency of the pulsar, its time deriva- 
tive of order i and the braking index. 



( V 0)) 2 



(1) 



We obtained n>20. Large braking indices between glitches 
are common among Vela-like pulsars and are likely to be as- 
sociated with glitch recoveries. These large values should not 
be interpreted as the long-term braking index due to secular 
spin evolution but instead these correspond to transient states 
caused by large gl itch activity as discussed in section 3.2.2 of 
Hobbs et al. (2010). Dipole braking indices have been measured 
onl y for a few pulsars wit h the highest spindown rates (see Table 
1 of lEspinoza et all d201 lh ). 

With its characteristic age of -21 kyr, PSR J 1856+0245 is a 
Vela-like pulsar possibly affected by large glitch activity, which 
we have not yet directly seen. Thus, we fix here the pulsar brak- 
ing index to the canonical value of 3. 

We assume the existence of three primary photon fields 
(CMBR, far IR (dust), and st arlight) and use the interstellar ra- 
diation field from GALPROP dPorter et al]|2 005) to estimate the 
photon fields at the Galactic radius of PSR J1856+0245. A dis- 
tance of 9 kpc in the direction of the pulsar corresponds to a 
Galactic radius of 5.4 kpc. At this radius, the peak of the SED 
of dust IR photons corresponds to a black-body temperature of 
T ~ 32 K with a density of ~ 1.1 eVcrrT 3 , while the SED of 



http://galprop.stanford.edu/resources.php?option=data 
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stellar photons peaks at T ~ 2500 K with a density of ~ 1.2 
eV cm 4 . 

A simple exponentially cutoff power-law injection of elec- 
trons, evolved properly over the pulsar lifetime, often provides 
an adequate match to PWN SEDs. For this injection spectrum 
we fit four variables: final magnetic field Bf, electron high en- 
ergy cutoff E cut , electron power-law index p, and initial pulsar 
spin period Pq. The best fit parameters, with errors, are given in 
Table 2. The initial spin period gives an age of 20 + 2 kyr. This 
model poorly matches the low energy MAGIC points, as shown 
in Figure [2] (Top). The low value of ~ 4 fiG found for the mag- 
netic field is consist ent with 6 nG found fo r the 15 kyr old PWN 
HESS J1640^I65 in lLemiere et all (l2009i) 

Another option to fit the multi-wavelength data is to adopt 
the relativistic Ma xwellian plus pow er-law tail electron spec- 
trum proposed by Smtkovskv (2008). We implement this spec- 
trum as described in Grondin et al.l(l201 lb . The best fit, presented 
in Fig. |2] (Middle) and Table 2, is obtained with kT = 0.64 TeV 
corresponding to an upstream Lorentz factor of 2.5 x 10 6 . A very 
high cutoff of 390 TeV is required, though the power-law index 
of p = 2.44 is co nsistent with the value of ~ 2.5 proposed by 
Spitkovskvl (120081) . The initial spin period gives an age of 15 + 2 
kyr. The relativistic Maxwellian plus power-law model better 
matches the multi-wavelength data, and also directly probes the 
upstream pulsar wind via fitting of the upstream Lorentz factor 
of the wind. 

A hadronic scenario is also possible, with y-rays arising 
from proton-proton interactions. For this model, corresponding 
to Fig. [2] (bottom) and Table 3, we fix the ambient gas density 
at 50 cirT 3 and age at 20 kyr. It should be noted though that the 
total energy injected in such scenario is very high even for the 
large density assumed here. 

6. Conclusions 

Using 3 years of Fermi-LAT data, a y-ray source has been de- 
tected at high significance at a position coincident with the TeV 
source HESS J 1857+026. We have investigated whether a model 
in which a PWN is powered by the pulsar PSR J1856+0245 
can reproduce the multi-wavelength data. In such a scenario, the 
VHE spectrum observed by MAGIC and H.E.S.S., combined 
with the limits imposed by the steep LAT data, is difficult to 
match with a simple power-law injection of electrons (or pro- 
tons), and we find a significantly better fit with a relativistic 
Maxwellian plus power-law spectrum. The low magnetic field 
of these leptonic fits, due to the stringent X-ray upper limit, im- 
plies that if PWN leptons are indeed responsible for the y-ray 
flux, they must be dominated by relic electrons which have es- 
caped the PWN core into weakly magnetized surroundings. The 
hadronic scenario relaxes this constraint, though the energy re- 
quirements are quite high even for a dense ambient medium, and 
a very hard power-law index is required. At present the true na- 
ture of HESS J1857+026 remains a mystery, though the new 
LAT data and X-ray upper limit hint that this source may be an- 
other relic PWN, increasing the population of such high energy 
y-ray systems. 
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Fig. 1. TS maps computed by pointlike. The green crosses 
show the sources of the 2FGL catalog included in the model, 
whereas the blue X represents the source we added to the 
model. The green contours show the H.E.S.S. data (Aharonian 
et al., 2008). The magenta circle gives the position of PSR 
J1856+0245. Top: TS map obtained between 10 and 300 GeV. 
The position of the Fermi-LAT excess is consistent with that 
of H.E.S.S. Note that HESS J1857+026 is not included in the 
model. Middle: TS map obtained between 0.1 and 1.3 GeV, 
showing the residual excess taken into account in our model. 
Bottom: Residual TS map obtained between 10 and 300 GeV 
when all sources are included. 
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Fig. 2. Spectral energy distribution of HESS J1857+026 with 
a simple exponentially cutoff power-law electron spectrum 
(Top), a relativistic Maxwellian plus power-law electron spec- 
trum (Middle), and an exponentially cutoff power-law pro- 
ton spectrum (Bottom). The X-ray flux upper limit obtained 
using Chandra (green), LAT spec tral points (red), MAGIC 
points (violet) dKlepser et alj|201lt) . and H.E.S.S. points (blue) 
(Ahar onian et al . 2008) are shown. The black line denotes the to- 
tal synchrotron, inverse Compton and pion decay emission from 
the nebula.Thin curves indicate the Compton components from 
scattering on the CMB (long-dashed), IR (medium-dashed), and 
stellar (dotted) photons. 
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Table 2. Model Parameters : the leptonic fits, including the statistical uncertainties on the fitted parameters. 



Model 


n" C v(kF) 


B(Bcc^) B f QiG) 


p(E-P) 


kT (TeV) 


E cut (TeV) 


P (ms) 


* 2 /d.o.f. 


Power-Law 
Rel. Max. 


3 C 3 C 
3 C 3 C 


-l.5 c 3.9 ±0.4 
-1.5 C 2.9 ±0.6 


2.12 ±0.03 
2.44 ±0.10 


0.64 ± 0.09 


120 ± 44 
390 ± 280 


11.1 ±6.8 
42 ±3 


21.9/21 
12.5/20 



" Pulsar braking index 

* Sedov phase onset, after which r oc r 03 

c Held fixed 



Table 3. Model Parameters : the hadronic fit, including the statistical uncertainties on the fitted parameters. 



Model rf t Sedm (kyr) 


B(Bcc^) Bf(nG) 


p(E-P) 


n(cm 3 f 


E cut (TeV) £ (xl0 50 erg) d 


* 2 /d.o.f. 


Hadron 3 e 3 e 


-1.5« 20™ 


1.83 ±0.04 


50 c - e 


75 ± 25 0.64 ± 0.06" 


24.6/21 



" Pulsar braking index 

b Sedov phase onset, after which r oc t 03 

c Ambient medium density 

d Total energy injected 

e Held fixed 



